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Abstract 

We update Monte Carlo simulations of the three-dimensional SU(3) + adjoint Higgs theory, 
by extrapolating carefully to the infinite volume and continuum limits, in order to estimate 
the contribution of the infrared modes to the pressure of hot QCD. The sum of infrared con- 
tributions beyond the known 4-loop order turns out to be a smooth function, of a reasonable 
magnitude and specific sign. Unfortunately, adding this function to the known 4-loop terms 
does not improve the match to four-dimensional lattice data, in spite of the fact that other 
quantities, such as correlation lengths, spatial string tension, or quark number susceptibilities, 
work well within the same setup. We outline possible ways to reduce the mismatch. 
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1. Introduction 



Motivated primarily by experimental heavy ion collision programs at RHIC, LHC and FAIR, 
there are on-going large-scale efforts to determine the pressure of hot Quantum Chromody- 
namics (QCD) through numerical Monte Carlo simulations of the full four-dimensional (4d) 
three-flavour [l]-[3] or even four-flavour [4J theory. Due to the significant cost of simulating 
light dynamical fermions near the chiral limit, these efforts typically make use of so-called 
staggered quarks. At any finite lattice spacing, the staggered action does not possess the 
flavour symmetries of the continuum theory, and is also problematic in the flavour-singlet 
sector. Although these problems will eventually be overcome by the use of less compromised 
fermion discretizations, it would appear welcome in the meanwhile to explore complementary 
avenues as well, in order to offer timely crosschecks on the results that are being produced. 

One potentially useful avenue in this respect is provided by effective field theory meth- 
ods. In the limit of a high temperature, T, and a small gauge coupling, g, QCD develops 
a hierarchy of three momentum scales, vrT ^ gT ^ g^T/tr. The largest scale, ttT, can be 
systematically integrated out, yielding a so-called Dimensionally Reduced effective field the- 
ory [5l[6], or "Electrostatic QCD" (EQCD) [7]. EQCD has previously been used to determine 
non-perturbatively quantities such as spatial correlation lengths [8], the spatial string ten- 
sion [9]-|llj. and quark number susceptibilities [T2j. In all cases, a surprisingly good match 
to the results of 4d lattice simulations was found, even down to temperatures very close to 
that of the deconfining crossover, (see, e.g., refs. [131 El)- We would therefore like to 
explore the extent to which a similar success could be achieved for the pressure (within pure 
Yang- Mills theory), and subsequently perhaps apply the same methods to physical QCD. 

A way to apply EQCD to the study of the pressure on the non-perturbative level, as well 
as first numerical results, were put forward a number of years ago [H]. The practical results 
of ref. [H] suffered from two problems, however: certain 4-loop logarithmic terms, whose 
form was not understood at the time, were effectively missed [16]; and the systematic error 
from the continuum extrapolation carried out was probably underestimated, given that the 
approach to the continuum limit has turned out to be more delicate than anticipated |17j . 
The purpose of this paper is, therefore, to update the analysis of ref. [15j, both by making 
use of the new theoretical ingredients in refs. [161 H^j. and by increasing the numerical effort 
manifold. Unfortunately, the match to 4d lattice data does not improve despite all these 
efforts. On the other hand, the fact that systematic errors are now under control, implies that 
the discrepancy needs to be taken seriously, and this offers us the possibility to speculate on 
the kind of physics that might be missing in our approach. In particular, we wish to discuss 
the role of higher dimensional operators within the EQCD framework, as well as the new 
qualitative features that should be expected from more radically improved effective theories. 

We start by specifying the general setup of our approach (Sec. [2]); go on to discuss the 
details of the lattice formulation within EQCD (Sec. [3]); present the main results (Sec. H]); 
and conclude with a discussion and outlook (Sec. [5]). 
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2. General setup 



At high temperatures, the pressure of QCD can be written as [7] 

PQCd(T) = Phard(r) +PsoftiT) ■ (2.1) 

Here Phard is a matching coefficient (defined in the MS scheme) which gets contributions only 
from the hard scale, k ~ ttT, and is computable in perturbation theory, while 

Psoft(r) = l^lim^ ^ In 1 VA'^VA^Q exp(-5E) |_ , (2.2) 

where V = Jd'^x. is the d-dimensional volume (d = 3 — 2e), represents the contributions of 
the soft scales. The effective action can be written as 

= Jd''x{^^-Tr[F^^]+Tv[D,,Aof + mlTv[Al] + X,{Tr[Al]f + ...Y (2.3) 

Here = {i/gs)[Di, Dj], Di = di - igsAi, Ai = AfT", Aq = A^T"", and are hermitean 
generators of SU(3). In the following we set e — > in all finite quantities, whereafter the 
dimensionalities of and A3 are GeV. 

Now, psoft is scale-dependent, just like Phard- Within the truncated form of Eq. (j2.3p . the 
scale dependence can be worked out explicitly [18]. Defining the dimensionless ratios 
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and ignoring terms of 0{g^) in terms of 4d power counting, we can write 
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where Jl is the MS scale parameter; is by definition the vacuum energy density of the 



(truncated) EQCD, computed with = g^ and scaled dimensionless by dividing with g^; and 
dA = iVc ~ 1) = -^c, -^c = 3. The function T'jj^ can, in turn, be written as 
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Here [5q = —0.2 it 0.8 is the non-perturbative "Linde term" from the pure three-dimensional 
Yang- Mills theory |19l I20j . estimated numerically in refs. [2T]-[23], while is the 4-loop 

perturbative contribution sensitive to the adjoint Higgs field [l8], specified for complete- 
ness in appendix A. The remainder, is what we address in the following. 

For dimensional reasons, the function necessarily depends on the parameter y (it gets 
contributions from five loops and beyond and, before the rescaling with g^, therefore contains 
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Figure 1: The stripes in the (x, ?/)-plane corresponding to various 4d theories in the temperature range 
T = (0.5 . . . 1500) A-j^ j24], with Ni denoting the number of massless quarks. The points in Table 1 
correspond to A'f = and Ni ~ 2. The relation to the temperature is y(jvj. = o) — 0.32 \log^Q{T / K-^) + 
0.91], y(^, . 2) ~ 0.38 [logio(r/A^) + 0.98], y^^^ = 3) - 0.39 [logio(r/A^) + 1.04]. 

terms of the form g^'^'^^ / m]^^ ^ with n > 0). Thus no information is lost (i.e. no integration 
constant is needed) if we take a partial derivative with respect to y, leading to a condensate: 

a,j^U-. y) = (T.' \mh - (TV - (TV ' P'**' 

where Aq = Aq/q^. For a non-perturbative study, we need to change the scheme from MS to 
the lattice, and then the relation in Eq. (12. Sp goes over into 

dyT^ix,y) = lim{(TV [ig]) - fo{m) - h{m) - y-\ ^(m) - y'^ h(jn)\ , (2.9) 

where 

and fi are functions that have recently been determined numerically in ref. [T7]. The task, 
therefore, is to measure (Tr [^oDa lattice, for values of x and y corresponding to the 

physical finite-temperature QCD (cf. Fig. [T]), and insert then the result in Eq. (12. 9p . 
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3. Lattice simulations 



The lattice study is carried out with the standard Wilson discretized action, 
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^Tr[io(x)t/.(x)io(x + i)^;(x)] 



+ E{«Tr[ig(x)] + ^l^(Tr[ig(x)] 



(3.1) 



where ?7j(x) is a link matrix; x + i 
and 



X + aii, with ii a unit vector; Pij(x) is the plaquette; 

[3^4-- (3-2) 
The bare mass parameter is given by [2_5j 
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(3.3) 



For historical reasons, we have implemented in most of the runs a partial 0{a) improve- 
ment [26], by taking 

0.328432 - 0.835282X + 1.167759a;2 



(3.4) 



but for the mass parameter no improvement was carried out, since the additive 0{a) terms 
have not been determined: yiatt = U- The improvement of xjatt turns out to play little 
practical role if properly taken into account in the subtraction of Eq. (|2.9p : in fact we did 
not implement it in the last sets of runs, corresponding to A'^f = 2 as well as to /? = 240 with 
A'f = 0. In any case our setup conforms with ref. [17|, where the functions fi in Eq. ()2.9p 
were determined numerically for general xiatt- 

The input parameters in any run are x, y, (3. The phase diagram of the system in this 
space has a "disordered", or symmetric phase, as well as a symmetry broken phase [27]. The 
simulations we carry out represent the physical QCD only in the symmetric phase [21]. It 
turns out that on the physical stripe (Fig. [1]) the symmetric phase is actually metastable, 
but strongly enough so for any practical effects to be miniscule in large enough volumes (see, 
however, the discussion in Sec. [5|). 

The three-dimensional SU(3) -|- adjoint Higgs theory is a confining gauge theory, and 
possesses a mass gap (for a practical demonstration see, e.g., ref. [3)- Thereby all finite 
volume effects must be exponentially suppressed. For completeness we have checked this 
explicitly: the condensate (Tr[AQ])a as a function of [3/N = where L is the extent 

of the box, is shown in Fig. [2j No finite- volume effects are visible in the range P/N < 1, or 
L > 6 /g^ , to which we restrict in the following. The parameters used for the production runs 
are listed in Table (H 
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Figure 2: Finite- volume values for (Tr [Aq])^, as a function of the physical extent P/N = G/g^L of the 
box {L = aN), for a small y (left) and large y (right). No volume dependence is visible for (3/N < 1.0. 
At small y and small volumes, the metastability of the physical phase is too weak to hold the system 
there for any length of Monte Carlo time. 



4. Numerical results 

In Fig. Owe show (Tr [^oDa values of y, as a function of the lattice spacing 1/(3 = ag'^/6, 

and after the subtraction of the various terms in Eq. (j2.9p . The figure indicates that after 
the subtractions, a continuum limit (1//3 0) can be taken; and that perturbation theory 
does converge even at the smallest y (corresponding to the lowest temperature), in the sense 
that each subtraction is smaller than the previous one, and that the remainder is at most of 
the same order as the last subtraction. Increasing y, perturbation theory converges faster, 
and the significance loss due to the subtractions becomes substantial. 

In Fig. [4] we show a magnification of the fully subtracted results, for selected values of y. 
This highlights the regular- looking approach to the continuum limit. 

In order to carry out the continuum extrapolation, we have tested four different fit func- 
tions: a constant fit, a fit linear in 1//3, a fit including a linear term and the logarithm 
(In /?)//?, as well as a quadratic fit. The results are shown in Fig. [5l We observe that all 
dependences apart from the first one lead to consistent results; the situation is thus quite 
different from ref. [12] , where fits allowing for a logarithmic term were the only ones leading 
to a sensible outcome. In other words, the 3- loop and 4- loop subtractions in Eq. (j2.9p . which 
were not available at the time of ref. [12], have effectively removed any possible logarithmic 
terms. In the following, intercepts obtained from linear extrapolations (squares) will be used. 
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Table 1: The continuum parameters x,y (cf. Eqs. ()2.4|) . ()2.5|) ): the lattice coupHngs /3 (cf. Eq. (j3.2|) ): 
and the box sizes N {V ^ a^N^) used for the production runs. Values marked with a star correspond 
to N{ = 2, others to Nf = 0. In the few cases where the box size has the superscript 2, two independent 
runs were launched. In total, our sample consists of 186 lattices. 

The dashed line in Fig. [5] shows the curve 

(Tr[^o]>i.-^(ci + ^ + f), (4.1) 

with ci = 0.0200(8), C2 = -0.0191(14), C3 = 0.0149(6). This representation describes our 
data well (x^/d.o.f. = 2.6) in the whole y-range considered. 

The next task is to integrate Eq. (|2.9p in order to determine the function J^^- Since 
Eq. (j2.9p contains a partial rather than a total derivative, and on the physical stripe the 
parameter x changes (cf. Fig.[T]), an integration is strictly speaking not possibl^; however, the 
leading x-dependent term of corresponds to a contribution of the type Psoft ~ 93-^3/''^3 ~ 

^In ref. [H] this problem was tackled by also measuring the condensate dxJ-j^ — {{Tr[Ag])'^)a — ... . 
However, some of the renormalization constants needed for the MS conversion remain unknown, and in any 
case the practical effects from this condensate appeared to be too small to change the qualitative behaviour. 
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Figure 3: Shown is (Tr [AQ])a at a small and large y, after the consecutive subtraction of the functions 
/o, /i, /2, fs of Eq. ()2.9p . corresponding to 1-loop, 2-loop, 3- loop and 4-loop effects, respectively. 



g^T^i which is of higher order than other effects we have ignored. More importantly, all effects 
containing the parameter x are numerically subdominant for Nc = 2> (cf. Eq. (lA.ip ). Thereby, 
close to the physical stripe, the remainder in Eq. ()4.ip can simply be integrated to 
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(4.2) 



We note that the numerical value of the integral is only weakly dependent on the ansatz 
used for the continuum extrapolation. 

We can now compare our result with 4d lattice data. In the following we consider pure 
glue [28], for which systematic errors are best under control. To get the EQCD result, we 
need to evaluate both phard and psoft (Eq. ()2.ip ). In practice, we take part of the terms in 
Psoft) namely those specified explicitly in Eqs. ()2.6p . ()2.7p . and combine them with the hard 
contribution Phardi specified in Eqs. (5), (6) of ref. [29j. Given that the non-perturbative con- 
stant contains numerical errors and that Phard contains a perturbative constant, denoted 
by Ahard in ref. [29], which remains unknown, we treat the combination A^ard + ^aC'I/Jg as 
a free parameter. The remaining soft contribution, 5psoft = —Tg^T^s (cf- Eq. (j2.6p ). is given 
by the sum of Eqs. ()4.2p . (jA.ip . The free parameter A^ard + dAC\(3G is fixed by minimizing 
the x^-difference of 4d lattice data [28] and our full result in the range T/Tc > 3.0. The 
outcome of the fit is shown in Fig. [6l (It corresponds to A^ard 

+ dAC\PG = -7.24.) 

It can be observed from Fig. [6] that our high-temperature result and the low-temperature 
4d lattice result depart already at T ~ 3.2Tc, and that the fit is in general not particularly 
good. In fact, compared with the fit in ref. [29j which assumed J-^ = 0, increases about 
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Figure 4: A magnification of the fully subtracted results from data of the type in Fig. [3l for selected 
values of y, indicated on the right. Multiple data points at the same /? correspond to independent 
runs, or different volumes (cf. Table [T]). 



30-fold. (The fit of ref. [29] corresponds to Ahard + ^aCI/Jg = -1-84.) The reason for the 
increase is that (Tr [Ag])^ in Eq. (j4.ip is positive, whereby in Eq. (j4.2ll is negative, 



whereby Psoit/T'^ gets a contribution increasing rapidly at small temperatures (cf. Eq. (j2.6p ). 
This is completely unlike the behaviour of the 4d lattice data. 

Now, one possible reason for the mismatch could be that, starting at 0{g'^), higher order 
operators should be added to the EQCD action [16j. It seems, however, that in practice such 
operators cannot change the result in a substantial way. Indeed, for A'^f = 0, the higher order 
operators were determined in ref. [30], and the only one contributing at 0{g'^T^) reads 
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Treating this perturbatively (as a "vertex"), and working in dimensional regularization, yields 
a contribution 

(IaCa g^ml 1 dAC\ 7 4 
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(47r)4T 47r 30 {Ai^f 

where we inserted Ca = 3 and ~ gT. In contrast, the effect that we have determined 
numerically in this paper has the magnitude (cf. Eqs. (j2.6p . (j4.2p ) 
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from fits with different /3-dependences, as indicated in the legend. Some points have been slightly 
displaced for better visibility. The dashed curve is from Eq. (14. ip . 

This is larger than Eq. (j4.4p by more than two orders of magnitude. Another effect of 0{g'^) 
is that the parameter y should be computed at the 3-loop order; however, this only modifies 
the relation of y and T/T^, and cannot revert the upward trend of our data at small T/T^. 
Therefore, it seems unlikely that infrared sensitive effects (from momenta k ~ gT) describable 
with simple improvements of EQCD would be the cause for the mismatch; the problems are 
perhaps more likely related to the treatment of the ultraviolet modes {k ~ vrT). 




5. Conclusions and outlook 

The purpose of this paper has been to approximate the non-perturbative contribution of 
the dynamics represented by a (truncated) effective theory called EQCD, Eq. (|2.3p . to the 
pressure of hot QCD. The result is constituted by the sum of the non-perturbative Linde 
term (Eq. ()2.7p ). perturbatively known terms up to 4-loop level (Eqs. p.6p . (jA.ip ). as well 
as an ultraviolet finite all-orders remainder that we have estimated numerically (Eq. (j4.2p ). 

On the technical side, the main content of our study was to carefully carry out the contin- 
uum extrapolation needed for estimating the remainder, Eq. (|4.2p . As has been illustrated 
in Fig. [5l the continuum extrapolation appears now to be under control, thanks partly to 
the recent determination of the functions fi in Eq. (|2.9p |17| . Therefore, as we have argued, 
the representation in Eq. (14. 2p should be free of substantial systematic errors close to the 
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Figure 6: Left: the best fit of our result to 4d lattice data for p/T^ at iVf = 28J. Our result comes 
out as a function of T/A^; for converting to T/Tc we scanned the interval Td Ky^ = 1.10 . . . 1.35 
(dark band). For comparison, with the light band we show the outcome for = Right: the 
corresponding trace anomaly, (e — 3p)/T^ — Td{p/T^)/dT. 



physical stripe of Fig. [TJ Comparing with the lower order terms in Eq. (jA.ip . Eq. (j4.2p also 
has a reasonable magnitude, and indicates that perturbation theory within EQCD is in fact a 
useful tool at all parameter values corresponding to the 4d theory. In principle, the numerical 
result could also be compared with improved resummation methods defined within EQCD 
(see, e.g., refs. [3T|-|33j). 

At the same time, from the phenomenological point of view, our result is somewhat of a 
disappointment: as shown in Fig. [6l the match to 4d lattice data in the range T/Tc > 3.0 is 
not particularly smooth. In fact, as shown in the figure, including the newly determined > 5- 
loop remainder decreases the quality of the fit significantly. This implies that, unfortunately, 
we are not in a position to realize our original goal of offering consolidated crosschecks for 
the A^f 7^ QCD pressure in the interesting temperature range T/Tc ~ 1-5 . . . 3.0. 

On the other hand, our study raises the theoretical question of what kind of effects could 
be responsible for the mismatch between our results and that of 4d lattice simulations. On 
the mundane side, one possibility would be a substantial contribution from the condensate 
denoted by dxJ~Ms = ((Tr [^o])^)a ~ ••• • Unfortunately its systematic inclusion would 
require a significant amount of new analytic and numerical work; moreover, as order-of- 
magnitude estimates and previous preliminary simulations suggest, it appears unlikely that 
this condensate could significantly change the qualitative behaviour that we have observed. 

On a more adventurous note, let us point out that, qualitatively, the reason for the mis- 
match is that the condensate in Eq. (j4.ip is too large, and grows rapidly as y (or T) de- 
creases (note that for A^f = 0, the range T/Tc = 10'^. ..10^ corresponds to y ~ 0.3. ..1.2, or 
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1/y ~ 3. 3. ..0.8, cf. Fig. [T]). Since the condensate measures the mean squared fluctuation of 
the field Aq^ this means that the colour-electric gauge field fluctuates too much. Indeed, if 
the remainder in Fig. [5] continued to decrease also for 1/y > 0.5 (in principle the remainder 
can even become negative), then the match in Fig. [6] could conceivably improve dramatically. 
Perhaps one reason for the large fluctuations at small y could be that the physical phase of 
EQCD is merely metastable [21] (see also Fig. 3 of ref. [27j)? If so, improved effective theories 
of the type in refs. j3l]-[39], which do not make any explicit distinction between the scales 
k ^ gT and k ^ ttT and which are by construction free of the metastability problem, might 
yield smaller fluctuations and a better outcome. 

Of course, if no scale separation is made between the scales k gT and k ~ vrT, then 
one should in principle also account for corrections of 0{g^) to the hard part of the pressure 
(Phard) Eq. (|2.ip ). rendering the analysis very hard. Nevertheless, it might be interesting to 
explore if a phenomenologically successful recipe could be found even without this last step. 
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Appendix A. Four- loop three-dimensional vacuum energy density 



The function !F^^"^, defined through Eq. ()2.7p and representing the contribution of the 
adjoint scalar Aq to the vacuum energy density of EQCD, with fl = 03 and with A'c = 3 so 
that only a single scalar self-coupling appears, can be written as |18pl : 



+ 



+ 



+ 



__dA_lP 
'(47r) 3 

dA y 



Ca 



(47r)2 4 

dA 
(47r)3 2 

W{dA + 2 

dA 



3-2 ln(4y) 
.2 



+ x{dA + 2)^ 



vr- 11 , 

— + In 2 

12 3 3 



+ - CA{dA + 2) 



21n(4y) - 1 



(4^) 



3 - ln(16y) 
43 491 



^^{dA + 2f] 



1536 



311 43 



vr^ ln(4y) + — + 



256 32 



ln2 



^We have inserted the value 710 = 7r^/24 — ln^2/2 [30], not known analyticaUy in ref. [18j . 



11 



+ llln^2 



461 2^ 491 2i ^ 1793 



9216 



1536 



5x _,9 



vr 



512 

17 



vr 



96 

43 11 



vr 



ln^(4y) -ln(4y) \ In 2 



X 



+— CA{dA + 2) 



6 



vr 



ln2(4y)+ (9- — )ln(4y) 



31 



8 + 61n2 + ^vr2 + ^ln2 + 81n2 2- ^C(3) 
24 4 4 ' 



21 



X 



CA{dA + If 



21n(4y) + 1 



ln(16y) - 1 



+^(dA + 2)(dA + - 



vr 



-vr^ ln(2y) + — - 21C(3) 



35 



1 ln(4y) + 1 + 3 In 2 vr^ + — In 2 C(3) 
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(A.l) 



References 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 

[7] 
[8] 
[9] 
[10] 

[11] 
[12] 

[13] 

[14] 



Y. Aoki, Z. Fodor, S.D. Katz and K.K. Szabo, JHEP 01 (2006) 089 |hep-lat/0510084l . 



C. Bernard et al, Phys. Rev. D 75 (2007) 094505 [ hep-lat/0611031] . 

M. Cheng et al, Phys. Rev. D 77 (2008) 014511 [0710.0354]. 

M. Cheng et al, PoS LAT2007 (2007) 173 [0710.4357]. 

P. Ginsparg, Nucl. Phys. B 170 (1980) 388. 

T. Appelquist and R.D. Pisarski, Phys. Rev. D 23 (1981) 2305. 

E. Braaten and A. Nieto, Phys. Rev. D 53 (1996) 3421 [hep-ph/9510408l . 

A. Hart, M. Laine and O. Phihpsen, Nucl. Phys. B 586 (2000) 443 [hep-ph /0004060|. 

B. Lucini and M. Teper, Phys. Rev. D 66 (2002) 097502 [hep-lat / 0206027] . 



P. Giovannangeh, Phys. Lett. B 585 (2004) 144 [hep-ph/0312307] 



M. Laine and Y. Schroder, JHEP 03 (2005) 067 (hep-ph/0503061]. 
A. Hietanen and K. Rummukainen, JHEP 04 (2008) 078 [0802.3979]. 
M. Cheng et al, Phys. Rev. D 78 (2008) 034506 [0806.3264]. 
Y. Maezawa et al, PoS LAT2008 (2008) 194 [0811.0426]. 



12 



[15] K. Kajantie et al, Phys. Rev. Lett. 86 (2001) 10 f hep-ph/0007109 |. 

[16] K. Kajantie et al, Phys. Rev. D 67 (2003) 105008 |hep-ph/021132l] . 

[17] F. Di Renzo et al, JHEP 09 (2008) 061 [0808.0557]. 

[18] K. Kajantie et al, JHEP 04 (2003) 036 [ hep-ph/0304048| . 

[19] A.D. Linde, Phys. Lett. B 96 (1980) 289. 

[20] D.J. Gross, R.D. Pisarski and L.G. Yaffe, Rev. Mod. Phys. 53 (1981) 43. 

[21] A. Hietanen et al, JHEP 01 (2005) 013 (hep-lat/04120l)8] . 

[22] A. Hietanen and A. Kurkela, JHEP 11 (2006) 060 [hep-lat/06090T5] . 

[23] F. Di Renzo et al, JHEP 07 (2006) 026 [h ep-ph/0605042l . 

[24] K. Kajantie et al, Nucl. Phys. B 503 (1997) 357 |hep-ph/9704416l . 

[25] M. Laine and A. Rajantie, Nucl. Phys. B 513 (1998) 471 [hep-lat/9705003] . 

[26] G.D. Moore, Nucl. Phys. B 523 (1998) 569 |hep-lat/970905l] . 

[27] K. Kajantie et al, JHEP 9811 (1998) Oil |hep-lat/9811004| . 

[28] G. Boyd et al, Nucl. Phys. B 469 (1996) 419 [hep-lat/9602"007] . 

[29] M. Laine and Y. Schroder, Phys. Rev. D 73 (2006) 085009 [ hep-ph/0603048] . 

[30] S. Chapman, Phys. Rev. D 50 (1994) 5308 [hep-ph/94073 13|. 

[31] J.P. Blaizot, E. lancu and A. Rebhan, Phys. Rev. D 68 (2003) 025011 [ hep-ph/0303045] . 

[32] M. Inui, A. Niegawa and H. Ozaki, Prog. Theor. Phys. 115 (2006) 411 ihep-ph/0501277] . 

[33] G. Cvetic and R. K5gerler, Phys. Rev. D 75 (2007) 054016 |h ep-ph/0612130] . 

[34] A. Vuorinen and L.G. Yaffe, Phys. Rev. D 74 (2006) 025011 [hep-ph/0604100] . 

[35] R.D. Pisarski, Phys. Rev. D 74 (2006) 121703 |hep-ph/060 8242]. 

[36] A. Kurkela, Phys. Rev. D 76 (2007) 094507 [0704.1416]. 

[37] A. Dumitru and D. Smith, Phys. Rev. D 77 (2008) 094022 [0711.0868]. 

[38] Ph. de Forcrand et al, Phys. Rev. D 77 (2008) 125014 [0801.1566]. 

[39] CP. Korthals Altes, 0810.3325 [hep-ph]. 

[40] Y. Schroder, unpublished; E. Bejdakic, 0810.3097 [hep-ph]. 



13 



